Thermodynamics of quantum gases for the entire range of temperature 
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We have analytically explored thermodynamics of free Bose and Fermi gases for the entire range 
of temperature, and have extended the same for harmonically trapped cases. We have obtained 
approximate chemical potentials of the quantum gases in closed forms of temperature so that the 
thermodynamic properties of the quantum gases become plausible specially in the intermediate 
regime between the classical and quantum limits. 
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I. INTRODUCTION 

In the undergraduate courses on statistical physics one 
has to plot specific heats of free (ideal homogeneous) 
Bose and Fermi gases Jj. To plot specific heat of a 
quantum gas with respect to its temperature, one has 
to know temperature dependence of its chemical poten- 
tial which in general, is obtained from an implicit rela- 
tionship with the number density of particles jl|, Q ■ This 
relationship involves a polylog function of an exponential 
of the chemical potential which can not be obtained in a 
closed form of temperature as because inverse of a poly- 
log function, in general, does not exist. In spite of that, 
students are taught how to get approximate temperature 
dependence of the chemical potential for the classical and 
quantum (degenerate) regimes But, for the interme- 
diate regime, they are advised to get graphical solutions. 
Feeling difficulties of obtaining the graphical solutions, 
we in this paper, will incorporate an easy technique of 
obtaining approximate chemical potentials of the quan- 
tum gases in closed forms of temperature not only for 
the intermediate regime but also for the entire range of 
temperature. 

After the observation of Bose-Einstein condensation in 
1995 {■} •") . studying trapped quantum gases now-a-days 
becomes a topic of high experimental and theoretical in- 
terest of research 6-8]. So, it is reasonable to extend 
our discussion for harmonically trapped (ideal inhomo- 
geneous) cases. 

Calculations of this article will begin with the Bose- 
Einstein statistics. In subsection-2A we will present a 
technique of how to get an approximate chemical poten- 
tial of a free Bose gas in a closed form of temperature. 
We will do the same for a harmonically trapped Bose gas 
in subsection- 2B, and for free and harmonically trapped 
Fermi gases in section-3. We will plot approximate chem- 
ical potentials in FIG 1. We will also plot specific heats, 
equation of states, and susceptibilities in FIGs 2 and 3 
using our approximate chemical potentials. 



II. THERMODYNAMICS OF BOSE GASES 

Let us consider a 3-D free Bose gas in thermodynamic 
equilibrium with a (heat and particle) reservoir which is 
characterized by temperature T and chemical potential 
fi. All the thermodynamic quantities of this system can 
be obtained from the Bose-Einstein statistics 
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where fii is the (equilibrium) average number of particles 
at the ith single particle state having energy e^. From 
Eqn.(l) we can see that total average number of parti- 
cles {N = ^ - fii) and total average energy {E = ^ ■ n^ei) 
are functions of fi. Hence, all the thermodynamic quanti- 
ties like pressure (p) , average number density of particles 
{n = N/V), average energy per particle (e), specific heat 
per particle (cy), etc. are expected to be functions of /i. 

For a free gas, we can replace the single particle energy 
eigenstates {i} by the single particle momenta {p} so 
that Ci essentially becomes Cp = /2m where m is the 
mass of a single particle. Degeneracy of the level Cp in 

the semiclassical limit is given by where V is 

the volume of the system. In the thermodynamic limit 
(A^ — > oo,y oo,N/V = const.), standard text book 
result for specific heat per particle of the free Bose gas is 
given by [l| 
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for T <T,. 
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is the Bose-Einstein conden- 



where = ^( 

sation temperature at and below which /i takes the high- 
est possible value (0), z — e''/'"'-^ is the fugacity, and 
\j\j(^z) — z + ~ + ~ + ... is a polylog function of order j. 
Plotting Cy for T < Tc is very easy. But, plotting c„ for 
T > Tc is not an easy job until one manages to get the 
temperature dependence of z or of ^ from the implicit 
relation [l| 
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which can be obtained as a result of phase-space integra- 
tion of the right side of Eqn.(l), and where t = T/Tc- 
As inverse of a polylog function (Lij(z)) does not exist 
except for j = 1, we are not able to get /i from Eqn.(3) 
or Cy from Eqn.(2) as a function of temperature without 
approximation in particular for T > Tc- 

On the other hand, for a harmonically trapped case, all 
the particles are 3-D harmonic oscillators, and the single 
particle energy levels are given by f-i — + i)hu}, where 
u! is the angular frequency of oscillations. Although the 
degeneracy (gi) of this level is i^/2 -I- 3i/2 -t- 1 (§], yet in 
the thermodynamic limit {N — >■ oo, w — & Nw^ = 
const.), only the first term of the degeneracy contributes 
significantly. Zero point energy can also be neglected in 
this limit. The implicit relation between the chemical 
potential and temperature (in the thermodynamic limit) 
for this case is given instead of Eqn.(3) by 0, i, ^ 
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where i = ^ and = ^[^J^^^ 

In the following we will present approximate analyses 
of Eqns.(3) and (4), and will obtain approximate chemi- 
cal potentials of the two cases indicating suitable suffices 
for different temperature regimes. 



A. For free Bose gas 

Classical regime of a quantum gas is characterized by 
t ^ 1 (or by z <^ 1) so that in this regime we may 
write Li3/2(2;) ~ z which along with Eqn.(3) yields the 



^3/2 

approximate classical result 



^kT-ln{t/CHm)- 



(5) 



On the other hand, the quantum regime for Bose gas, 
is characterized by t < 1 (or by z ^ 1). For t < 1, we 
have ^ = 00, and for t > 1, ^ is to be obtained from 
Eqn.(3). In the following we will extend a standard (but 
not much familiar) result [ll| for t ^ 1 up to the second 
lowest order. From Eqns.(l) and (3) we can redefine the 
polylog function Li3/2(z) as 
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in an integral form where v ~ \fi\/kT. One can check 
that Li3/2(2:) is not an analytic function of z particularly 
about z — 1. Thus to know even the approximate values 
of chemical potential for t ^ 1 is a tricky job. Let us 
adopt the trick of Landau and Lifshitz for this purpose 



[13 . Subtracting U^/M) 
the above equation as [12| 
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from 'Li^/2{z) we can recast 



e^de. (7) 



1 



It is clear from Eqn.(7) that, for t ^ 1 (or for v ^ 1), 
the integrand will contribute only for smaller values of e, 
so that up to the second lowest order we can have 
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Let us use this result at the left hand side of Eqn.(3). In 
its right hand side, we may have the Taylor expansion 
about t = 1 as 

lA=l-(i-l) + (<-l)2 + .... (9) 
Plugging above two expansions in to Eqn.(3) we get 
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up to the second lowest order uy 5t — t — I. It should 
be mentioned that the lowest order term in Eqn.(lO) was 
previously obtained by Bhattacharjee adopting the same 
trick of Landau and Lifshitz [ll|, [l^l ■ 

Intermediate regime is an overlap of the classical and 
quantum regimes, and it can be characterized by i ^ 1. 
For this case, our approximation is made considering only 
the first three terms of Li3/2(-2) instead of its all terms in 
Eqn.(3). Thus solving the cubic equation of z we have 
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where /i(t) = [[36\/6 - <dV2]t^/^ + lUt^CiD + 
\/3([3072\/3 - 864]i9 + [3456^6 - 864V2]C(|)t^^/2 + 

mi2t^C{l)Y'^Y'^ . Eqn.(ll) eventually is not only ap- 
plicable for the intermediate regime but also for the clas- 
sical regime. 

For t < 1, we have iJL[t) = 0, and have nothing to 
be approximated. For t > 1, we have two independent 
approximate results in Eqns.(lO) and (11). These two 
results intersect at t k, 1.711. Thus for i > 1, we get 
the approximate chemical potential as a combination of 
these two in the following form 
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where is a unit step function. Now it is time to verify 
how good is our approximate chemical potential. We can 
verify it comparing with the (exact) graphical solutions 
of the chemical potential. 

Let us now outline how to get the graphical solutions. 
In a figure we plot left hand side of Eqn.(3) with respect 
to z, and right hand side of Eqn.(3) for t — ti,t2,t^ etc. 
Intersecting points of these plots are the solutions {{zi}) 
of z for t = ti,t2,t3 etc. Thus we can have a number 
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(a): For free Bose gas (b): For trapped Bose gas (a): For free Bose gas (b): For trapped Bose gas 

p/kTj /j/kT, c,,/k c,,/k 




FIG. 1: Solid lines in FIGs la, lb, Ic & Id represent ap- 
proximate chemical potentials in Eqns.(12), (16), (23) & (26) 
respectively. Dotted and dashed lines represent classical and 
quantum results respectively. Points represent exact graphi- 
cal solutions of the respective chemical potentials. 



of graphical solutions {fii — kTiln{zi)} for given values 
of t, and plot them in FIG. la. In the same figure we 
also plot our approximate chemical potential in Eqn.(12), 
and see that it matches reasonably well with the exact 
graphical solutions. Hence, we can approximately plot 
other thermodynamic variables using our approximate 
chemical potential in Eqn.(12). As an example, we plot 
specific heat per particle (from Eqn.(2)) in FIG. 2a, and 
compare with its exact graphical solutions. 

That the approximate Cy matches well with the exact 
graphical solutions gives us courage to adopt our approx- 
imation technique to apply for the harmonically trapped 
ideal Bose gas, free Fermi gas, and for the harmonically 
trapped ideal Fermi gas as well. 



B. For trapped Bose gas 

Chemical potential of the harmonically trapped Bose 
gas is to be obtained from Eqn.(4). The classical re- 
sult (/i>(i) ~ — 3/cTln[t/C3 (3)]) for this case can be im- 
proved for the intermediate regime considering only the 
first three terms of Li3(z) in Eqn.(4), and it results 



9 leSi'^ 



(13) 



18 



where /2(t) = [261t^ + 256i6C(3) + 16\/2(458i 
261ii5C(3) + 128^12(2 (3)^ i/Y/3_ the other hand, in 
the quantum regime in particular for t ^ 1, we should 
expand Li3(z) and 1/t in Eqn.(4) about z = 1 and t = 1 
respectively. While the expansion of Li3(z) can be ob- 
tained (using Li3(z) — ^(3) = '^'^J^-' dx and expanding 



FIG. 2: Solid lines represent approximate specific heats per 
particle. Points represent their respective exact graphical so- 
lutions. Dotted lines represent classical results. 



Li2(x) in a manner of obtaining Eqn.(8)) as 



Li3(z)-C(3)-C(2)z^ 



C(2)^2 
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that of 1/t has already been obtained in Eqn.(9). Pres- 
ence of the logarithmic term in above expansion says that 
Li3(2:) is nonanalytic about z = 1. However, plugging the 
expansions in Eqns.(14) and (9) in to Eqn.(4) we get [l3| 



-3^(i-l) 



(15) 



up to the lowest order in t — 1. Once again for t > 1, we 
have two independent approximate results in Eqns.(15) 
and (13). These two results intersect at t ~ 1.065. Com- 
bining these two we get our desired result for < > 1 as 



^>(t) w ^>(t)6'(1.065 -t) 



1.065). (16) 



We plot this approximate chemical potential in FIG. lb, 
and compare it with its exact graphical solutions. 

In the following section we will extend our discussions 
for the Fermi systems. Although in the classical regime, 
behavior of a Fermi gas is identical with that of a Bose 
gas, yet in the quantum regime, they are quite differ- 
ent. This difference needs special care for extending our 
approximation technique to be applied for the Fermi sys- 
tems. 



III. THERMODYNAMICS OF FERMI GASES 

For an ideal Fermi gas, average number of particles 
occupying ith single particle state is given by the Fermi- 
Dirac statistics 
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(17) 
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While for a Bose gas we have defined t as T/Tc, hence- 
forth for a Fermi gas, let us redefine t to be T/Tp where 
Tp is the Fermi temperature. Integrations of right hand 
side of Eqn.(17) result [T] 

(-Li3/2(-^)r(5/2))'/' = l (18) 
for a free Fermi gas with Tp = ^llfe (S"'^"')^^'^ [IIj and 

(-Li3(-z)r(4))'/' = i (19) 
for a harmonically trapped Fermi gas with Tp — 

In the following we will present approximate analyses 
of Eqns.(18) and (19), and will obtain approximate chem- 
ical potentials of the two cases indicating suitable suffices 
for different temperature regimes. 

A. For free Fermi gas 



(a): For free Bose gas (b): For free Fermi gas 
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(c): Pauli paramagnetism (d): Landau diamagnetism 

A^Paiili^f/C ^Landau ^f/C 




FIG. 3; Solid lines in FIG 3a and 3b represent equation of 
state (p) in units of p''^'' = fikT. Solid lines in FIG 3c and 
3d represent paramagnetic (xPauii) and diamagnetic (xLandau) 
parts of susceptibilities of the ideal homogeneous Fermi (elec- 
tron) gas in units of C/Tp = fJ-onfi%/kTF. Points represent 
respective exact graphical solutions. 



Adopting the previous manner we get the classi- 
cal result for the free Fermi gas to be /i>(i) ~ 
— |A;Tln[tr3 (5/2)], and improved it considering only the 
first three terms of — Li3/2(— z) to get the approximate 
result for the intermediate regime {t > 1) as 
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where fsit) = [l927rt3 
\/3(l22887r2t6 -|- [1152\/2 - 

3072V3]7r3t9)^^Y/^ On the other hand, for the quan- 
tum regime (t < 1) as well as for z ^ 1, we ex- 
pand — Li3/2(— -z) in Eqn.(18) according to Sommerfeld's 
asymptotic formula 1] 
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n=l ^ 

{j-[2n-l]) 
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Plugging the above expansion for j 
we get [1J| 



/i<(<) = kTp 
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B. For trapped Fermi gas 

Adopting similar treatment obtaining approximate re- 
sults of the free Fermi gas, we get the classical and im- 
proved classical result for the trapped Fermi gas respec- 
^2Qt^vely as M»(t) w -3fcrin[tr^(4)] and [l| 

"9 165 X 3i/3/4(i) 32/3 - 



(21) 



3/2 in to Eqn.(18) 



(22) 



For the entire range of temperature, we have two in- 
dependent approximate results in Eqns.(22) and (20). 
These two results intersect at f « 0.723. Their combi- 
nation in the following form 

fi(t) w fi<{t)e{0.72S -t)+ fi>it)e{t - 0.723) (23) 

gives us our desired approximate chemical potential. We 
plot this approximate chemical potential in FIG. Ic, and 
compare it with its exact graphical solutions. 
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8/4 (t) 



(24) 



where Uit) = t/[128 - 783^^ + 16(64 - 783^^ + 
8244t^)^/2]^/'^. For the quantum reg ime, using the 
asymptotic formula for j = 3 we get [l3| 



fi<{t) = kTp 



(25) 



Once again, for the entire range of temperature, we have 
two independent approximate results in Eqns.(25) and 
(24). These two results intersect at t « 0.468. Their 
combination in the following form [l5| 

ti{t) « ^<(i)6i(0.468 -t) + fi>{t)0{t - 0.468) (26) 

gives us our desired approximate chemical potential. We 
plot this approximate chemical potential in FIG. Id and 
compare it with its exact graphical solutions. 



IV. CONCLUSIONS 

Although it is impossible to get chemical potentials of 
the quantum gases in exact closed forms of temperature, 
yet we have obtained them in approximate closed forms. 
The closed forms give plausibility to reproduce thermo- 
dynamics of the quantum gases for the entire range of 
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temperature. As a claim of plausibility we have plotted 
specific heats of all the quantum gases in FIG 2, and 
equations of states (pressures) of the free quantum gases 
in FIGs 3a & 3b, and paramagnetic & diamagnetic parts 
of the susceptibilities of the ideal homogeneous Fermi gas 
of spin 1/2 in FIGs 3c & 3d using the respective approx- 
imate chemical potentials in FIG 1. 

We have not shown relationships between specific heats 
and fugacities for the plots in FIGs 2b, 2c and 2d. These 
are easily derivable and available in Refs.[l3], [H and [ll] 
respectively. We have not also shown the derivations of 
equations of states and that of the susceptibilities in FIG 
3. These are available in standard text books p], fllj. 

For harmonically trapped systems, we have considered 
oscillations to be isotropic. Otherwise angular frequency 
(a;) in our results has be replaced by the geometric mean 
of clIj;, i^y, and ujz 2, 6]. 

It should have been mentioned that an approximate 
form of specific heat of the free Bose gas was previously 
obtained for T ^ Tc by Wang using Robinson's method 
which involves Mellin transform as well as analytic con- 
tinuation of Lij(z) [l^ [13], and for T ^ Tc by London 
using virial expansion method [11, [l^. These methods 
(or techniques) can also be used for the harmonically 
trapped Bose gas. On the other hand, approximate form 
of thermodynamic variables of the free and harmonically 
trapped Fermi gases can be obtained for T <Tp using 
Sommerfeld's asymptotic expansion of — Lij(— -z), and for 
T ^Tp using the virial expansion method [l|, [^ . While 
none of the above methods/techniques covers the entire 
range of temperature, our approximation technique cov- 



ers the entire range not only for the Bose gases but also 
for the Fermi gases. However, it would be an interest- 
ing problem for undergraduate and graduate students to 
get a temperature dependent form of specific heat of the 
trapped Bose gas for T in the way we have got ap- 
proximate chemical potentials using the trick of Landau 
and Lifshitz [l^. Students can also generalize the virial 
expansion to get approximate temperature dependence 
of specific heats of the trapped Bose and Fermi gases. 

To conclude, this paper illustrates approximate anal- 
yses of polylog functions in demonstrating temperature 
dependence of thermodynamic quantities of ideal quan- 
tum gases for the entire range of temperature in an easy 
and plausible way. In this paper we have applied our 
approximation technique only to a few thermodynamic 
properties of the ideal quantum gases. Our technique 
can also be applied for other thermodynamic properties 
(eg. temperature dependence of entropy, energy, etc.) of 
the ideal and weakly interacting quantum gases [T5| . All 
the calculations of this paper have been done within the 
scope of undergraduate and graduate students. They can 
extend our approximation technique to get better results 
including higher order corrections, and generalize our ap- 
proximation technique to other spatial dimensions. 
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